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Abstract 

We consider a new generalization of the skew-normal distribution introduced by 
Azzalini (1985). We denote this distribution Beta skew-normal (BSN) since it is a 
special case of the Beta generated distribution (Jones (2004)). Some properties of 
the BSN are studied. We pay attention to some generalizations of the skew-normal 
distribution (Bahrami et al. (2009), Sharafi and Behboodian (2008), Yadegari et al. 
(2008)) and to their relations with the BSN. 

1 Introduction 

The skew-normal distribution (SN), introduced by Azzalini (1985), has been studied and 
generalized extensively. The aim of this article is to introduce a new family of distributions, 
which generalizes the skew normal, that is flexible enough to support both unimodal and 
bimodal shape. This new family, called Beta skew- normal (BSN), arises naturally when we 
consider the distributions of order statistics of the skew-normal distribution. The BSN can 
also be obtained as a special case of the Beta generated distribution (Jones (2004)). In 
this work we pay attention to three other generalizations of the skew-normal distribution: 
the Balakrishnan skew-normal (SNB) (Balakrishnan (2002), as a discussant of Arnold and 



Beaver (2002), Gupta and Gupta (2004), Sharafi and Behboodian (2008)), the generalized 
Balakrishnan skew-normal (GBSN) (Yadegari et al. (2008)) and a two parameter general- 
ization of the Balakrishnan skew- normal (TBSN) (Bahrami et al. (2009)). The above three 
extensions are related to the Beta skew-normal distribution for particular values of the pa- 
rameters. 

Given a random sample Xi, - ■ ■ , X„ from a distribution F{x), in general the distribution of 
the related order statistics does not belong to the family of F{x). In this paper we show 
that the maximum between the X^s from a Balakrishnan skew-normal with parameters m 
and 1, denoted by Xj ~ SNBm{l), is still a Balakrishnan skew-normal with parameters k 
and 1, where /c is a function of m and n. 

This paper is organized as follows: after describing briefly, in Section 2, the skew-normal 
distribution, its generalizations and listing their most important properties, in Section 3 we 
present some generahzations of the Beta distribution. In the last Section we define the Beta 
skew-normal distribution, we present its properties and some special cases. In particular 
the BSN contains the Beta half-normal distribution (Pescim et al. (2010)) as limiting case. 
Besides, we investigate its shape properties. We derive its moment generating function and 
we also compute numerically the first moment, the variance, the skewness and the kurtosis. 
We present two different methods which allow to simulate a BSN distribution. We explore 
its relationships with the other generalizations of the skew-normal and we show that the 
distributions of order statistics from the skew-normal distribution are Beta skew-normal and 
are log-concave. Furthermore, in this Section we give some results concerning the SNB dis- 
tribution. In particular we derive the exact distributions of the largest order statistic from 
SNBjn{i) and the shortest order statistic from SNBjn{—l). 

2 The skew-normal density and its generalizations 

The present Section recalls some important definitions and properties about the skew- 
normal distribution and some of its extensions. 
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2.1 The skew- normal density 

The skew-normal distribution refers to a parametric class of probability distributions 
which includes the standard normal as a special case. A random variable Z is said to be 
skew-normal with parameter A, if its density function is 

(^{z]\) = 1(^{z)^{\z) with A,2eM (1) 

where and $(■) denote the standard normal density and distribution, respectively. We 
denote a random variable Z with the above density by Z ~ SN[X). The parameter A controls 
skewness. The standard normal distribution is a skew- normal distribution with A = 0. We 
remind some properties of the SN distribution. 
Properties of SN(\): 

a. As A — )■ oo, A) tends to the half-normal density. 

b. If Z is a SN{X) random variable, then —Z is a S'A^(— A) random variable. 

c. If Z ~ SN{X), then Z^ ~ 

d. The density ([T]) is strongly unimodal, i.e. log 0(2;; A) is a concave function of z. 
The corresponding distribution function is 

/z pXt 
/ <p{t)<p{u)dudt = $(2) - 2T{z, A), 
-00 J —00 

where T(z, A) is Owen's function. The properties of this function are: 

1. -T(z,A)=T(z,-A); 

2. T{-z,\)=T{z,\)- 

3. 2T{z,l) = ^{z)(!>{-z). 

Using the properties of Owen's function, we have immediately the following ones: 
Property 1.1- $(-2;; A) = ^{z; -A). 
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Property 2. ^{z;l) = ^{z^. 

Property 3. $(2; A) + $(2; -A) = 2^{z). 

The class of skew-normal distributions can be generalized by the inclusion of the location 
and scale parameters which we identify as ^ and ip > 0. Thus if X ~ SN{X) then Y = ^ +ipX 
is a skew-normal variable with parameters ^, ip, A. We denote F by F ~ SN{C,, ip, A). 

2.2 The Balakrishnan skew-normal density and its generalization 

Balakrishnan (2002) proposed a generalization of standard skew-normal distribution as 
follows: 

Definition 1. A random variable X has Balakrishnan skew-normal distribution, denoted by 
SNBn{X), if it has the following density function, with G N, 

fn{x; A) = c„(A)0(a;)$(Aa;)", x eR, A G M. (2) 

The coefficient c„(A), which is a function of n and the parameter X, is given by 

Cn(A) 



where U ~ iV(0, 1). 

For n = and n = 1, the above density reduces to the standard normal and the skew- 
normal distribution, respectively. 

For n = 2 a random variable X with the density ^ is denoted by X ~ NSN(X) with 
C2(A) = arctan (Sharafi and Behboodian (2006)). 

The class of Balakrishnan skew-normal can be generalized by the inclusion of the location 
and scale parameters which we identify as /i and cr > 0. Thus if X ~ SNBn{X) then 
Y = fi + crX is a Balakrishnan skew-normal variable with parameters /i, a, X. We denote Y 
by Y SNB„ifi,a,X). 
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Remark 1. Sharafi and Behhoodian (2008) have shown that for A = 1, (|2| is the density 
function of the {n + l) — th order statistic X(^n+i) a sample of size n+1 from iV(0, 1). 
Moreover, for A = — 1, ^ is the density function of the first order statistic in a sample 
of size n+1 from A^(0, 1). 

Recently, Yadegari et al. (2008) introduced the following generalization of the Balakrish- 
nan skew-normal distribution and explained some important properties of this distribution. 

Definition 2. A random variable X is said to have a generalized Balakrishnan skew-normal 
distribution, denoted by GBSNn,m{^) , with parameters n, m G N and A G M, its density 
function has the following form: 

fn,U^; A) = —^<P{xMXxr (1 - $(Aa;))'" x G M, (3) 

where C„,^(A) = EI^o (?) (-1)7-!! 

For m = this density reduces to the Balakrishnan skew-normal. 

Remark 2. Let Xi, ■ ■ ■ , X„ be a random sample from a N{0, 1) . Then the j-th order statistic 
is a G B S N j_i.a-j{l) , with j = 1, ■ ■ ■ ,n. In this case we have that 

/ -x poo I 

c,-^,n-,{i) = E 7 ' J (-1)' y_ <i>ixmxy-'^^dx = ^^_,^l^_^y - (4) 

Bahrami et al. (2009) discussed a two-parameter generalized skew-normal distribution 
which includes the skew-normal, the Balakrishnan skew-normal and the generalized Balakr- 
ishnan skew-normal as special cases. 

Definition 3. A random variable Z has a two-parameter Balakrishnan skew-normal distri- 
bution with parameters Ai, A2 G M, denoted by T„^m(Ai, A2), if its pdf is 

fn,m{z- Al, A2) = j^—MzMXlZrHX2Zr, Z G M, (5) 

andn, m are non-negative integer numbers. The coefficient c„^m(Ai, A2), which is a function 
of the parameters n, m, Ai and X2, is given by 

Cn,Ml^ A2) = E [<l>(AiX)"<l>(A2X)-] , 

where X ~ A^(0,1). 
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The following properties are direct consequence of definition (|3] 
Properties of T„^m(Ai, A2): 



1. TBSNi,i{\i,0) = SN{\i) and TS^A^i,i(0, A2) = SN{\2); 

2. TBSNn,m{X,X) = SNBn+m{Xy, 

3. TBSNn,U^,,0) = SNB^iX,) and TS5iV„,„ (0, A2) = ^iV5„(A2); 

4. TBSN^^mi^ -A) = GBSN^^mW and TS^iV„,^(-A, A) = Gfi5iV^,„(A); 

5. TBSNn,UO, 0) = TBSNo,o{Xu A2) = iV(0, 1). 



Remark 3. Let Zi, - ■ ■ , Zn i.i.d N{0, 1) and < Z(2) < ■ ■ ■ < &e the corresponding 
order statistics, then Z(^i) ~ TBSNi_i^n-ii^, — !)■ 

The location-scale two-parameter Balakrishnan skew-normal distribution is defined as 
the distribution of F = /i + aX, where X ~ TBSNn,m{Xi, A2)- Hence /i G M and a > 0, are 
the location and scale parameters, respectively. We denote F by F ~ TBSNn,m{fJ', o", Ai, A2). 

In the rest of the paper we denote by A) the density function of SN{X) and by 
fn,m{-] Ai, A2) the density function of TBSNn^mi^i, ^2)- 

3 Some extensions of the Beta distribution 

In this Section we introduce two families of distributions which generalize the Beta one. 

3.1 The generalized beta of the first type 

We recall the following definition due to McDonald (1984). 

Definition 4. A variable X is said to have a generalized beta distribution of the first kind 
with positive parameters a,b,p and q if its density is given by 
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If p = 1 and q = 1 the variable X is a Beta of the first kind with parameters a and b. 
For q = 1 and a = 1 the variable X is said to have a Kumaraswamy distribution with 
parameters p and 6 (Kumaraswamy (1980)). 

3.2 The Beta generated distribution 

Starting from the beta distribution, Jones (2004) defined a new family of probability 
distributions, called Beta generated distribution. Following the notation of Jones, the class 
of beta-generated distributions is defined as follows. 

Definition 5. Let F{-) be a continuous distribution function with density function /(■)■ The 
univariate family of distributions generated by F{-), called beta generated distribution, with 
parameters a,b > 0, has pdf 



gF{x;a,b) 



1 



iFix)r-\l-Fix)f-'fix) 



(7) 



B{a,b) 



where B{a,b) is the complete beta function. 



Thus, this family of distributions has distribution function given by: 




a, b> 




where the function Ip^x) denotes the incomplete beta ratio defined by 




(9) 



with 




< y < 1 



(10) 



the incomplete beta function. Replacing ^ and (10) in ([s]), we get that this family of 



distributions has distribution function 
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Remark 4. Let /(■) unimodal and continuously differentiable, if a,b> 1 then gp{-,a,b) is 
also unimodal. The strong unimodality, i.e. log- concavity, of f{-) implies strong unimodality 
of gf{-,a,h). 

Eugene et al. (2002) studied in details the family of beta-normal distribution (BN) and 
discussed its properties. Recently Pescim et al. (2010) proposed the beta half-normal (BHN) 
distribution to extend the half-normal (HN) distribution. We now recall the definitions of 
Beta-normal distribution and Beta half-normal distribution. 

3.2.1 Beta-normal distribution 

When in ([T]) F{x) is the normal distribution function with parameters n and a we have 
the Beta-normal family with distribution function 

Gf(^)(x;a,6) = -^^^*^ ^ - zf-' dz (12) 

and corresponding probability density function 




where a~^(f) (^^) and $ (^^) are the normal density and distribution with parameters /x 
and a. 

A random variable X with Beta-normal distribution with vector of parameters ^ = (0, 1, a, b) 
is denoted by X ~ BN{a, b). 

3.2.2 Beta half-normal distribution 

Let F{x) = 2$(x) — 1, with x > 0, the distribution function of the half normal distribu- 
tion. By using F(x) in ([T]) the density function of the Beta half-normal distribution (BHN) 
is given by 

^2V)_i(x;a,6) = ^I— (2$(x) - 1)"-^(1 - ^x)f-'<j){x), x > (14) 
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and the relative distribution function is 

Gf<,(.)_i(x;a,6)= / ^^(2$(t)-ir-i(l-$(t)f-V(t)rft, a; > 0. (15) 

^ ' Jo B{a,b} 

When X is a random variable following the BHN distribution, it is denoted by X ~ 

BHN{a,b). 

4 A new generalization of the skew-normal distribu- 
tion: the Beta skew-normal 

In this Section we define the Beta skew-normal class and we present some of its proper- 
ties. Replacing in ([T]) F{x) by $(a;;A), we obtain the Beta skew-normal distribution, with 
distribution function given by 

Gf(,^,)(x; A, a, b) = z'^-\l - zf-^dz (16) 

and probability density function 

^7f(,.^,)(a;; A, a, b) = ^^^{Hx; X)r-\l - <!>{x; X)f-'<P{xMXx). (17) 

Throughout this paper, we denote the Beta skew-normal distribution with vector of param- 
eters ^ = (A, a, b) by BSN{X, a, b). The class of the Beta skew-normal can be generalized by 
the inclusion of the location and scale parameters which we identify as fi and a > 0. Thus 
if X ~ BSN{X, a, b) then Y = fi + aX is a Beta skew-normal with vector of parameters 
^ = (/i, cr. A, a, b). We denote Y by F ~ BSN{fi, a, A, a, b). 

We now present some properties concerning the BSN{X, a, b). 
Properties of BSN{X,a,b): 

a- 9iix;\)i^'^ ^> 1' 1) = ^) for all X eR, i.e. BSN{X, 1, 1) = SN{X). 

b. 5(|(^..o)(x;0,a,6) = g^^.^^{x;a,b) for all x eR, i.e. BSN{0,a,b) = BN{a,b). 

c. ^|(^.o)(a;; 0, 1, 1) = (p{x) for all x G M, i.e. BSN{0, 1, 1) = N{0, 1). 
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d. 1, 1, 1) = 0(x) for all x G M, i.e. BSN{1, 1) = iV(0, 1). 

e- -1' 1' = '/'l^^) for all x G M, i.e. 55iV(-l, 1, ^ = N{0, 1). 

f. If X ~ BSN{X, a, 6) then -X ~ 5^iV(-A, 6, a). 

g. If X ~ BSN{X, a, b) then F = $(X; A) is a Beta{a, b). 

h. If X ~ BSN{X, a, b) then F = 1 - <I)(X; A) is a 5eta(6, a). 

i. As A — )■ +00, fi'f (^.;^) (a;; A, a, b) tends to the Beta half-normal density. 

Remark 5. Properties from a to e establish that the family of BSN{X, a, b) contains the 
standard normal distribution, the skew-normal distribution and the Beta-normal distribution 
as special cases. 



Proof. The proof follows directly from (17) and from elementary properties of the skew- 



normal distribution. □ 

The BSN distribution is easily simulated using Property g as follows: if Y has a beta distribu- 
tion with parameters a and b, then the variable X = ^^^(Y; A) has BSN{X, a, b) distribution, 
where $~^(.; A) is the quantile function of the skew-normal distribution. In Figure 1 are plot- 
ted random samples generated by the BSN distribution for some a, b, X with the respective 
curve of the density function obtained using R. 
By Property f we can deduce the following proposition. 

Proposition 1. Let X ~ BSN{X, a,b) and Y ~ BSN{—X,b,a). We have the following 
statements: 

1. Ex{X) = -Ey{Y); 

2. varx{X) = varyiY); 

3. 71 (X) = -71 (y); 
I 72(X) = 72(1^); 
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with 7i and 72 we indicate the skewness and the kurtosis, respectively. 



Now we find the moment generating function of X wliich has density (17) 



Property 4. The moment generating function of X BSN{X,a,b) is given by 

Mxit) = -^^^e^Ez A))'^-^(l - \)f-'^{\Z)\ [li 

where Z ~ N{t, 1). 

We have the following recursion formula: 
Property 5. Let k and k>2. If X BSN{X, a, b), with a,b> 1 then 



Ex{X^) = {k- l)Ex{X'^-^) + XEx 



+ 



$(AX) 

+ (a + 6 - l)Eu [U^-^(t)(U; A)] - (a + 6 - l)Ev [V^-^(t){y] A 
where U ~ BSN{X, a — l,b) and V ~ BSN{X, a,b — 1) are independent random variables. 



Proof. The proof follows easily from application of the formula for integration by parts and 

dx 



by using the well note = —x(p(x) (see Arnold et al.(1992)). □ 



The Beta skew- normal density is in general asymmetric (see Figures 2 and 3). We have 
a partial result concerning symmetry: 

Proposition 2. If a = b and BSN{X,a,b) is symmetric about then A = 0. 
Proof. We consider the density of a random variable X ~ BSN{X, a, a): 

9iix;x)i-^; A, a, a) = ^-^0(x)$(-Aa;)(l - ^x; -A))'^-^(<l>(a;; -A))'^"^ 
^ ' ' B{a,a) 



this is equal to g^^^.^-j{x; X,a,a) if <I>(Ax) = $(— Ax) and ^{x;X) = $(a;; — A). However for 
Property [3] we find that $(a;; A) = ^{x) which implies that A = 0. □ 



Remark 6. Eugene et al. (2002) have shown that the BN{a, b) = BSN{0, a, b) is symmetric 
about when a = b. 
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Remark 7. From Remark^we know that, ifa,b> 1, the density (17) is strongly unimodal, 
i.e \ogg^f^^.y^-^{x\ A, a, h) is a concave function of x. We don't have general results for a and/or 
h < 1. A numerical study has shown that, when at least one of the two parameters a and h 
is closed to zero (0.10,0.20), the density can he bimodal (see figure 3). Numerically Eugene 
et al. (2002) observed that the BN is bimodal when both parameters a and b are less then 
0.214. 

Moments of the BSN cannot be evaluated exactly. We have computed them numerically 
using the software R. In Table (1) we have reported the values of the mean (ibsn, standard 
deviation cbsn, skewness 71 and kurtosis 72 for different values of the parameters a, b and 
A. From these numerical study we have noted that: 

• for fixed values of a and b the mean (ibsn and skewness 71 are increasing function of 
A; 

• for fixed values of b and A the mean ^bsn and skewness 71 are increasing function of 
a; 

• for fixed values of a and A the mean fiBSN is a decreasing function of b. 

We now give some results concerning the distribution of order statistics from a skew-normal 
distribution: 

Proposition 3. Let Xi, ■ ■ ■ ,X„ be a random sample from a SN{X). Then the j-th order 
statistic is a BSN{X,j, n — j + 1), where j = 1 ■ ■ ■ ,n. 

Proof. The proof readily follows using the standard formula of the density of X(j), the i-th 
order statistic of a random sample of size n from the distribution SN{X). □ 

From Proposition [3] follows immediately that the family of BSN contains the distributions 
of the order statistics of the skew-normal distribution. 
In particular we have the following corollaries: 
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CoroUciry 1. Let X 



,Xn be a random sample from a SN{1). Then 



= max{Xi, ■ ■ ■ 



is a BSN{l,n,l). 



Corollary 2. Let X 



,Xn be a random sample from a SN{—1). Then 



= mm{Xi, ■ ■ ■ ,Xn} 



is a BSN{-l,l,n). 

CoroUciry 3. Let < X(2) < • • • < be the order statistic from a sample of size n 
from a SN{\) distribution. Then -'^(i), i — 1, - • ■ ,n, has log-concave density. 

Proof. From Property d of Section 2 we know that Xi has a log-concave density. We conclude 
the proof using the following result due to Gupta (2004): Suppose < X(2) < • • • < 
be the order statistic from a sample of size n from a distribution having a log-concave density 
function. Then -^(i), i = 1, ■ ■ ■ ,n, has log-concave density. □ 

We now derive other properties of the BSN distribution. 

Theorem 1. Let X ~ BSN{\, a, b) be independent of a random sample (Yi, - ■ ■ , F„) from 



SN{X), then X\ {Y^„) < X) - BSN{X,a + n,b) and X\ (Y(i) > X) ~ BSN{X,a,b + n), 



where Y(„) and Y(i) are the largest and the smallest order statistics, respectively. 




(19) 



Also 



P(Y(„) <X)^P{Yi<X,--- ,Yn< X) 




B{a + n,b) 



B{a,b) 
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Taking derivative from (19) with respect to w, we obtain the BSN{X,a + n,b) density 
function, the proof of X\ (Y(i) > X) ~ BSN{\, a,b + n) is similar. □ 

Theorem 2. IfX ~ BSN{X, a, h) is independent of the random sample {Ui, ■ ■ ■ , Un-i, Vi, - ■ ■ , Vm-i] 
from SN{X) then X\ (t/(„-i) < X, V^i) > X) ~ BSN{X, a + n - l,b + m - 1), where U^n-i) 
and V(i) are the largest and the smallest order statistics, respectively. 

Proof. The proof is quite similar to the one of Theorem [T] □ 

Theorem [l] can be used to generate X ~ BSN(\, n, 1) by extending the acceptance- 
rejection technique, due to Azzalini, as follows (see Azzalini (1985) and Sharafi and Be- 
hboodian (2008)): first we generate a random sample T, f/i, f/2, ■ ■ ■ ,Un^i from SN{X), if 
max([/i, U2, ■ ■ ■ , Un-i) < T we put X = T. Otherwise, we generate a new random sample, 
until the above inequality is satisfied. 



4.1 Further results concerning the Balakrishnan skew- normal and 
the Beta skew-normal 

In this Section we present some results concerning the SNB distribution and link the 
distributions introduced in Section 2.2 with the Beta skew-normal. First we consider two 
results about the Balakrishnan skew-normal. We study the distribution of the largest order 
statistic from SNBm{l) and subsequently the distribution of the smallest order statistic from 
SNBm{—^)- We found that these distributions belong to the family of SNB. 

Proposition 4. Let Xi, ■ ■ ■ ,X„ be a random sample from a SNBm{l). Then 

X(„) = max{Xi, ■ ■ ■ ,Xn} 
is a SNBk{l), where k = n{m + 1) — 1. 

Proof. The proof follows easily using the standard formula for the density of X(„), the largest 
order statistic of a random sample of size n from the distribution SNBm{l)- □ 

In particular the following corollaries hold: 
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Corollary 4. Let Xi, ■ ■ ■ ,X„ be a random sample from a SN{1). Then 

= max{Xi, ■ ■ ■ ,Xn} 

IS aSNB2n^i{l). 

Proof. The skew-normal distribution with parameter A = 1 is a Balakrishnan skew-normal 
with parameters A = 1 and m = 1. 

The same result can be established making use of the well-known result for the density of 
the largest order statistic from the distribution SN{1) and Property [2] □ 

Corollary 5. Let Xi, ■ ■ ■ ,X„ be a random sample from a SNBm{—l)- Then 

= min{Xi, ■ ■ ■ ,X„} 

is a SNBk{—l), where k = n{m + 1) — 1. 

It follows immediately from corollaries ([T]), ([2]), (|4]) and ^ that the BSN distribution 
is related to the previous skew-normal generalizations. In fact, its density simplifies to the 
Balakrishnan skew-normal when 6 = 1, A = 1 and a > 1 integer (or a = 1, A = — 1 and 
6 > 1 integer). Further, if A = it reduces to the generalized Balakrishnan skew-normal 
when a and b are both integers. These consideration have been summarized in the following 
proposition. 

Proposition 5. The BSN distribution satisfies the following properties: 

• = f2n-iAx] 1,0) for all xeR, i.e. BSN{l,n,l) = TBSN^n-iA^^O); 

• fi'f{x;-i)(^^ 1'"^) = fn,2m-i{x; 0, -1) for all X e R, I.e. BSN{-1, l,m) = TBSNn,2m-i{0, -1); 

• 5'|(a.;o)(^'0'"''"^) = 1,-1) for all X E R, i.e. BSN{0,n,m) = T5S'A^„_i_^„i(l, -1); 

where n, m are positive integer numbers. 

Given a random variable X ~ BSN{X, a, b) we are interested in constructing a random 
variable Y with Kumaraswamy distribution. This goal can be achieved using the below 
properties which follow easily from Property g and Property h of Section 4, respectively. 
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_l_ 

Property 6. If X ~ BSN{X,l,b) then Y = ($(X;A))" is a Kumaraswamy{a,b). In 
particular, if X ~ SNB2b-i{—l) then Y = {1 — $(— X)^)" is a Kumaraswamy{a,b). 

1 

Property 7. If X BSN{X,a,l) then Y = {1 — $(X;A))'' is a Kumaraswamy{b,a). In 
particular, if X SNB2a-i{l) then Y = {1 — is a Kumaraswamy{b, a) . 

Recently, Ferreira and Steel (2006) have presented a general approach which allows to 
generate skew distributions. They show that every univariate continuous skew distribution 
can be obtained from a "perturbation" of a symmetric one as it explained in the following 
definition: 

Definition 6. A distribution S is said to be a skewed version of the symmetric distribution 
F, generated by the skewing mechanism P, if its pdf is of the form 

s{y\f,p) = f{y)p[F{y)], yeR (20) 

where f and F are the pdf and cdf of a symmetric distribution on the real line, respectively, 
and p is the pdf of a distribution on (0, 1). 

Note that, if F is the standard normal distribution and p on (0, 1) is given by 

p{u; A, a, b) = -^A_$(A$-i(«)) ($($-i(«); A))"^"' (l - ^^-^u); A))'"' , (21) 



formula (20) reduces to Beta skew-normal with parameters X,a,b. Then the pdf of a Beta 
skew-normal with parameters A, a, b can be seen as a weighted version of (f){y), with weight 
function given by p ($(j/); A, a, b). 



16 



BIBLIOGRAPHY 

Arellano- Valle, R. B., and Azzalini, A. (2006). On the unification of families of skew-normal 
distributions. Scand. J. Statist, 33, 561-574. 

Arnold, B. C, Balakrishnan, N., Nagaraya, H. N. (1992). A first course in Order statistics. 
New York: Wiley. 

Arnold, B. C, and Beaver, R. J. (2002). Skewed multivariate models related to hidden 
truncation and/or selective reporting (with discussion). Test, 11(1), 7-54. 

Azzahni, A. (1985). A Class of Distributions Which Includes the Normal Ones. Scandinavian 
Journal of Statistics, 12(2), 171-178. 

Bahrami, W., Agahi, H., Rangin, H. (2009). A Two-parameter Balakrishnan Skew-normal 
Distribution. J. Statist. Res. Iran, 6, 231-242. 

Balakrishnan, N. (2002). Discussion of "Skewed multivariate models related to hidden trun- 
cation and/or selective reporting". Test, 11, 37-39. 

Cooray, K. and Ananda, M. M. A. (2008). A generalization of the half-normal distribution 
whit applications to lifetime data. Communication in Statistics & Theory and Methods, 37, 
1323-1337. 

David, H. A., Nagaraya, H. N. (2003). Order statistics. Hoboken, NJ: Wiley. 

Eugene, N., Lee, C, Famoye, F. (2002). Beta-normal distribution and its applications. 
Communications in Statistics - Theory methods, 31(4), 497-512. 

Ferreira, J. T. A. S, Steel, M. F. J. (2006). A constructive representation of univariate 
skewed distributions. Journal of the American Statistical Association, 101 (474) , 823-829. 

Gupta, R. C., Gupta, R. D. (2004). Generahzed skew-normal model. Test , 13(2), 501-524. 

Jamalizadeh, A., Balakrishnan, N. (2009). Order statistics from trivariate normal and 



17 



t,^— distributions in terms of generalized skew-normal and skew-tj, distributions. Journal 
of Statistical Planning and Inference, 139, 3799-3819. 

Jones, M. C. (2004). Families of distributions of order statistics. Test, 13(1), 1-43. 

Kumaraswamy, P. (1980). A generalized probability density function for double-bounded 
random processes. Journal of Hydrology, 46, 79-88. 

Loperfido, N. (2006). Statistical implications of selectively reported inferential results. 
Statistics & Probability Letters, 56, 13-22. 

Loperfido, N. (2008). A note on skew-elliptical distributions and linear functions of order 
statistics. Statistics & Probability Letters, 78, 3184-3186. 

McDonald, J. B. (1984). Some generalized functions for the size distribution of income. 
Econometrica, 52, 647-664. 

Pescim, R. R. Demetrio, C. G. B., Cordeiro G. M. E., Ortega, M. M.,Urbanoa, M. R., (2010). 
The beta generalized half-normal distribution. Computational Statistics and Data Analysis, 
54(4), 945-957. 

Sharafi, M., Behboodian, J. (2006). A new skew-normal density. J. Statist. Res. Iran, 3, 
47-61. 

Sharafi, M., Behboodian, J. (2008). The Balakrishnan skew-normal density. Statistical 
Papers, 49, 769-778. 

Yadegari, I., Gerami, A., Khaledi, M.J. (2008). A generalization of the Balakrishnan skew- 
normal distribution. Statistics and Probabilty Letters, 78, 1165-1167. 



18 



(a) 



(b) 



a 


6 


A 


fJ-BSN 




71 


72 


0.25 


0.25 


-10 


-1.1579 


1.4029 


-1.1329 


3.7648 






-1 


-0.6501 


1.9679 


-0.2378 


2.7777 












2.3382 


-0.0004 


2.6217 






1 


0.6484 


1.9649 


0.2306 


2.1362 






10 


1.1580 


1.4027 


1.1329 


3.7632 


0.25 


0.5 


-10 


-1.5906 


1.3469 


-0.7185 


2.7580 






-1 


-1.4424 


1.6716 


-0.3284 


3.0202 









-0.9631 


1.9061 


-0.0849 


2.8029 






1 


-0.1772 


1.5265 


0.0938 


2.8543 






10 


0.5446 


0.8728 


1.5054 


5.1988 


0.5 


0.25 


-10 


-0.5447 


0.8727 


-1.5061 


5.2003 






-1 


0.1773 


1.5265 


-0.0938 


2.8541 









0.9625 


1.9051 


0.0819 


2.7927 






1 


1.4411 


1.6694 


0.3203 


2.9849 






10 


1.6339 


1.3974 


0.8434 


3.2655 


0.5 


0.5 


-10 


-0.8979 


0.8874 


-0.9703 


3.3176 






-1 


-0.5882 


1.2659 


-0.1811 


2.9514 












1.5253 





2.8615 






1 


0.5882 


1.2659 


0.1811 


2.9514 






10 


0.9179 


0.9153 


1.0703 


3.7747 


0.5 


1 


-10 


-1.3018 


0.9148 


-0.8262 


3.4815 






-1 


-1.1664 


1.0704 


-0.3085 


3.1159 









-0.7043 


1.2479 


-0.1372 


2.9831 






1 





0.9999 





2.9999 






10 


0.4873 


0.5778 


1.3199 


4.8561 



a 


b 


A 


Hbsn 


cbsn 


71 


72 


0.5 


10 


-10 


-2.3678 


0.7314 


-0.7505 


3.7967 






-1 


-2.3617 


0.7389 


-0.7188 


3.7849 









-2.0809 


0.8033 


-0.6173 


3.5736 






1 


-1.0893 


0.6117 


-0.5642 


3.4799 






10 


-0.0182 


0.1429 


0.3706 


3.8635 


1 


0.5 


-10 


-0.4873 


0.5777 


-1.3200 


4.8570 






-1 





1 





3 









0.7043 


1.2479 


0.1372 


2.9831 






1 


1.1664 


1.0704 


0.3086 


3.1161 






10 


1.3018 


0.9148 


0.8262 


3.4814 


1 


1 


-10 


-0.7939 


0.6080 


-0.9556 


3.8232 






-1 


-0.5642 


0.8256 


-0.1369 


3.0617 












1 





3 






1 


0.5642 


0.8256 


0.1369 


3.0617 






10 


0.7939 


0.6080 


0.9556 


3.8232 


10 


1 


-10 


-0.0839 


0.1364 


-0.7082 


4.2018 






-1 


0.6744 


0.4536 


0.3597 


3.2722 









1.5388 


0.5868 


0.4099 


3.3314 






1 


1.8675 


0.5251 


0.5005 


3.4685 






10 


1.8807 


0.5124 


0.5744 


3.5243 


1 


10 


-10 


-1.8807 


0.5124 


-0.5744 


3.5243 






-1 


-1.8675 


0.5251 


-0.5005 


3.4685 









-1.5388 


0.5868 


-0.4099 


3.3314 






1 


-0.6744 


0.4536 


-0.3597 


3.2722 






10 


0.0839 


0.1364 


0.7082 


4.2018 



Table 1: The first moment, the standard deviation, the skewness and the kurtosis of 
BSN{X, a, b) for different values of a, b and A 
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;ure 1: Generated samples of the BSN distribution for some values of A, a and b 
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Figure 2: The Beta skew-normal BSN{X, a, h) for values of a,b > 1 
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Figure 3: The Beta skew-normal BSN{X, a, b) for values of a,b < 1 
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